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f"^ | Abstract 



We study the dynamic response of a [111] quantum impurity, such as lithium or cyanide 
in alkali halides, with respect to an external field coupling to the clastic quadrupolc 
moment. Because of the particular level structure of a eight-state system on a cubic 
site, the elastic response function shows a biexponential relaxation feature and a van 
Vleck type contribution with a resonance frequency that is twice the tunnel frequency 
A/%. This basically differs from the dielectric response that does not show relaxation. 
Moreover, we show that the elastic response of a [111] impurity cannot be reduced to 
that of a two-level system. In the experimental part, we report on recent sound velocity 
and internal friction measurements on KC1 doped with cyanide at various concentrations. 
At low doping (45 ppm) we find the dynamics of a single [111] impurity, whereas at 
higher concentrations (4700 ppm) the elastic response rather indicates strongly correlated 
defects. Our theoretical model provides a good description of the temperature dependence 
of Sv/v and Q^ 1 at low doping, in particular the relaxation peaks, the absolute values of 
the amplitude, and the resonant contributions. From our fits we obtain the value of the 
elastic deformation potential jt = 0.192 eV. 



1 Introduction 

The low-temperature properties of alkali halides may be significantly modified by the presence 
of substitutional impurities, such as Li, CN, or OH. Already a few ppm of these defects 
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totally change the thermal behavior and the elastic and dielectric response below Helium 
temperature. 

Regarding the polar molecules CN and OH, the point symmetry of the impurity site gives 
rise to several equivalent orientations; corresponding off-center positions arise for the small 
lithium ion. Quantum tunneling between these states results in a ground state splitting of 
about 1 Kelvin. Since the number of such tunneling states exceeds, even at low concentration, 
the number of small-frequency phonon modes of the host crystal, the impurities govern the 
low-temperature properties of the material. The cubic symmetry of fee crystals favors 8 defect 
positions in [111] -directions (CN and Li) or 6 positions in [100]-directions (OH in KC1). The 
resulting energy spectra have been discussed in detail by Gomez et al Q, and agree well 
with the Schottky peak observed by Pohl and co-workers for various impurity systems at low 
doping @j. 

In this paper we are concerned with the elastic susceptibiltiy of a [111] -impurity. Because 
of the cubic symmetry, such a tunneling system behaves in many respects as an ensemble 
of three two-level systems with energy splitting A || ||]. For example, the specific heat 
contribution of a [lll]-impurity is three times the Schottky peak of a two-state system. A 
similar relation holds true for the dielectric response function, since the dipolar transitions 
occur between adjacent levels only, as indicated by the dashed lines in the actual energy 
spectrum shown in Fig. Qj. This analogy has been used to describe lithium or cyanide defects 
in terms of a two-state approximation j|, ||, ||, [7| || . 

The elastic response function, however, shows a more complex behavior, resulting from the 
tensor character of the quadrupole operator Qij and the structure of the energy spectrum 
shown in Fig. |l[ According to the allowed quadrupolar transitions indicated by solid arows 
in Fig. |l|, sound waves or external stress perturb the impurity in two ways: First, they mix 
states separated by twice the splitting, giving rise to a van Vleck type susceptibility with 
resonance frequency 2A/h. Second, they lift the degeneracy of the states of a given triplet 
level, resulting in a relaxation contribution to the response function. 

Early measurements in KC1:CN samples || in the classical regime (k^T S> A) showed that 
the change of sound velocity of a T2 9 -mode (e = [110], k = [001]) varies as expected like 
\/T with temperature. Yet in addition, these experiments showed a change of sound velocity 
for an E^-mode (e = [110], k = [1-10]) with an amplitude of about ten times smaller than 
the T2 9 ones, whereas this mode should be unaffected according to the simplest model of a 
symmetric [lll]-tunneling defect. 

When measuring the sound velocity ( T2 9 -mode ) of KCl:Li as a function of temperature 



in the tunneling regime, Hiibner et al. found a hump at k-^T ~ A [10]. They related this 
observation to relaxational transitions between degenerate states, by calculating the static 
elastic response as the second derivative of the free energy of the impurity ion. Their results 
are valid in the low- frequency limit, where the external frequency is smaller than the impurity 
relaxation rates. 

The internal friction, however, cannot be obtained from a static theory, and the low-frequency 
limit is not always justified for the sound velocity. Starting from standard dynamic perturba- 
tion theory, we develop a dynamic theory for the elastic response function of a [lll]-impurity 
coupled to a phonon heat bath. The main purpose of this paper is to point out the pecu- 
liarities of the elastic response function that arise from the degenerate states of the energy 
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Figure 1: Energy spectrum of a [111] -impurity with zero asymmetry energy. Dashed arrows 
indicate the allowed dipolar transitions, and full arrows the quadrupolar ones. 

spectrum shown in Fig. |]. Moreover, we report sound velocity and internal friction data 
of a T2 S -mode for KC1 doped with cyanide molecules, and discuss them in terms of our 
[111] -impurity model with damping. 

The outline of our paper is as follows. In Sec. II we present the model and the basic equations 
for the elastic response function of a single impurity. The temperature dependence of the 
internal friction and sound velocity is evaluated in Sec. III. In Sec. IV we give some details 
of the experiments on KC1:CN, which are discussed in Sec. V. The final section contains a 
summary. 



2 Theory 

2.1 Quantum operators of a [lll]-impurity 

We start by discussing the properties of the energy eigenstates and by presenting our notation 
based on pseudospin operators. As shown in Fig. ||, the energy spectrum contains four 
equidistant levels; the upper and lower ones, Ai g and A2 U , are single quantum states, whereas 
the middle levels Ti u and T2 g are threefold degenerate. This groundstate multiplet can be 
described as the direct product of three two-level systems with energy splitting A. 

Each energy eigenstate is a superposition of states |r) localized at the eight impurity positions 



r = -d{al,a 2 z ,al) 



(1) 



with the effective two-state coordinates a\. For a lithium impurity, the off-center positions 
r form a cube of side-length d, whereas for cyanide impurities r indicates the orientation of 
the cigar-shaped polar molecule. Due to quantum tunneling along the axes i = 1,2,3, the 



energy eigenstates factorize as |cr^)|cr^)|cr^), where the variables a x = ±1 label odd and even 
superpositions \a l x ) = 2~ 1 ' 2 { \ri = d/2) + a x \ri = —d/2) } along the i-axis. 

Accordingly, the quantum operators of the impurity may be written as direct products of 
pseudospin operators a l a , where a = x,y,z label the usual Pauli matrices and a = the 
identity operator. In terms of quantum states at ri = ±d/2 we have 

a x = \d/2)(-d/2\+\-d/2)(d/2\ (2) 

al = \d/2)(d/2\-\-d/2)(-d/2\ (3) 

Thus a lithium or cyanide impurity is described as the product of three two-state variables. 
We adopt the shorthand notation for its quantum operators 

A aPl = a l a ®a}®a^ (4) 

where i = 1,2,3 label the crystal axes and Greek indices a = 0,x,y,z the pseudospin oper- 
ators. When denoting by ^A the tunneling amplitude between adjacent impurity positions, 
we obtain the Hamiltonian 

A 
H s = -— (A x qq + A 0xQ + A 00x ) . (5) 

Its eigenstates may be labelled according to the eigenvalues ±1 of a x ; for example, the 

ground state Ai g reads as |++- 1-), and the states of the triplet level Ti u as | — h+), \-\ h), 

and |+- 1 — ). According to the above discussion, the statistical operator p = e - ^ s /tr(e - ^ s ) 
factorizes as 

p = p 1 ®p 2 ®p 3 , p i = -(l + aiUnh(A/2k B T)y (6) 

(In the limit of zero temperature, the even state a x = 1 is occupied with probability 1.) From 
the statistical operator it is evident that the specific heat anomaly due to N such defects is 
identical to the Schottky peak of 3iV two- level systems with the Hamiltonian ^Aa x . As to 
the dynamic properties, a similar relation holds true for the dielectric susceptibility. This 
is not surprising, since dipolar transitions involve a single coordinate Vi and occur between 
adjacent levels with splitting A only; cf. Fig. ||. 

A different situation, however, is encountered when considering the elastic response: The 
quadrupole operator 

Qij =r i r j {l -6tj) (7) 

involves two coordinates, and hence causes transitions between energy eigenstates that differ 
in two labels a x . Thus elastic perturbations induce two types of transitions as indicated by 
the solid arrows in Fig. [l|. 

2.2 Phonon coupling 

Noting (|lj) and A z qqAq z q = A zz q, we may write the quadrupole operators as Q12 = (d/2) 2 A zz o, 
etc. When absorbing the factor (d/2) 2 in the coupling constant 7, we obtain the impurity- 
phonon interaction || 

Hsb = 2 2_^ 7 Q (A zz $ e" 2 + A z $ z e" 3 + A 0zz 623) • (8) 



When evaluating the strain tensor at the impurity site R and taking the limit of long wave- 
lenths, kd «C 1, we find 



Here k and e a are wave and polarisation vectors, and a labels the longitudinal and transverse 
phonon branches. 

2.3 Asymmetry energy 

Up to now we have considered an impurity in a potential that reflects the cubic site sym- 
metry. The strong dipolar and elastic interactions of nearby impurities, however, break that 
symmetry and give rise to an effective asymmetry energy, 

V = -^- A z00 - ^- A 0z0 - ^- A 00z , (10) 

where the Vi are random quantities with zero mean. We assume a Gaussian distribution with 
width a: 

p(v) = n p *(o with p ^) = 7^ exp (-\ S) (11) 

In principle, the Vi are dynamic variables closely related to the position operators of nearby 
impurities. In the case of weak interaction, Vi <C A, however, they may be treated as static 



random fields. Eq. (10) accounts for random fields arising from dipolar interactions; elastic 



coupling would lead to a potential that involves quadrupole operators such as A zz q. 

3 Elastic response function 

The lattice vibrations of the host crystal act on the tunneling impurity as a heat bath at tem- 
perature T. Upon a perturbation by an external strain field, the impurity regains the thermal 
equilibrium through absorption or emission of resonant phonons. Since the impurity-phonon 
coupling is weak, the resulting damping rates can be evaluated in second-order perturbation 
theory. 

3.1 Dynamic perturbation theory 

The linear response of the impurity to a time-dependent external strain field is described by 
the commutator correlation function {[Qij(t),Qij])- For convenience, we use dimenensionless 
quadrupole operators, e.g. A zz q = (2/d)Qi2 with (A 2 zzQ ) = 1; all relevant response functions 
can be expressed in terms of 

X (t) = ([A zz0 (t),A zz0 }) (12) 

It turns out convenient to calculate first the correlation matrix 

C a /3 Tt §et(t) = (A a /3ry(t)\Ag e £), (13) 



where we have defined the symmetrized correlation function 

(A\B)= l -{AB + BA). (14) 

Our evaluation of the correlation matrix is based on the Mori-Zwanzig projection formalism; 
the resulting memory kernel is calculated to second order in the impurity-phonon coupling 
[11]. Thus we obtain a matrix equation for the Laplace transform of (13), 



"W = ,_ n --' *(,) "■ (15) 

with the usual definitions of the static correlations 

M mn = {A m \A n ) = (r? _1 ) mri , (16) 



the frequency matrix 
and the memory kernel 



^^mn — \A rn \L,ApjTlp n , V^'J 



K mn {z) = {A m \(z - C ) ^Apjrjpn. (18) 

For convenience, we have replaced the triple index pair of @ with m, n, and we have used 
the shorthand notations A{ = i[HgB,Ai] and TiCq* = [Hs + Hb, *}, where Hb = J2 ^kb[pk- 
The thermal average (...) is with respect to p and the equilibrium distribution of the phonons. 

Solving the impurity dynamics now amounts to calculating the eigenvalues und residues of 
the resolvent [z — $7 — K\~ l . Since there are 4 3 operators A a ^, the matrix equation involves 
a 64-dimensional space. The matrices C, M, £1, and K being block-diagonal, the actual 
problem simplifies significantly. 

For symmetry reasons, the time evolution of the three operators A zz q, A z q z , and Aq zz is the 
same. Therefore it is sufficient to evaluate the dynamics in the invariant subspace containing 
one of them, e.g., 
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(19) 
In this invariant subspace, O, M, and K are represented by 8 x 8-matrices that are easily 



calculated according to (16-18). Note that the commutation relations for the operators Ai can 
be traced back to those for Pauli matrices. It can be shown that O and M can be diagonalized 
simultaneously; in the present example it is straightforward to find the corresponding unitary 
transformation U and to obtain 

O = U f nU = diag{0, 0, 0, 0, -2A, -2A, 2A, 2A)/R (20) 

M = U'MU = diag(m±, mi, rri2, 1^12, rri3,rn4, 1713,1714) (21) 



with 



mi = (1 - t) (1 - t 2 ) 
m2 = (l+t)(l-t 2 ) 
m 3 = (l +t ){l+t 2 ) 
m 4 = (l-t){l+t 2 ) 



where we have used j3 = l/k^T and t = tanh(/3A/2). 

Besides the fourfold degenerate zero frequency, there are two doubly degenerate finite fre- 
quencies ±2A/ft in Liouville space. While and M could be diagonalized simultaneously, 
this is not possible for the memory matrix K. It turns out, however, that its transform K 
is diagonal in each of the three degenerate subspaces of fi, with bare frequencies 0, ±2A/ft; 
the finite off-diagonal entries connect, e.g., the subspace of zero frequency eigenvalue with 
that of eigenvalue 2A/h. Because of K <C A/ft, we may neglect these off-diagonal entries of 
K; they would result in very small corrections to the frequencies and residues of the order of 
(ft-fC/A) 2 . In physical terms, the off-diagonal parts of K in the degenerate subspaces vanish, 
since they involve the phonon density of states at zero frequency; in a three-dimensional body, 
however, the density of phonon modes vanishes in the limit ujk — ► 0. 

Applying the usual Markov approximation, we evaluate the memory functions K(z) at the 
corresponding bare frequencies z = 0, ±2A/ft. It turns out convenient to separate real and 
imaginary parts according to K = K 1 + \K" . Whereas the entries in the subspace of zero 
eigenvalue are purely imaginary, those belonging to the finite frequencies ±2A/ft are complex 
numbers. Their real parts K'(±2A/h) would induce a shift of the resonance frequencies from 
the bare values ±2A/ft to slightly smaller ones. This shift being very small, we discard it 
and retain the imaginary, or dissipative, part K" only. The latter form a diagonal matrix 
f = QK"(z ) with z = 0, ±2A/ft. 

With these approximations for the memory kernel, £1, M, and T are diagonal; the correlation 
matrix thus factorizes, 

Cn(z) = ~ M] \~ , (22) 

Z 1 Ijj + 11 jj 

where j = 1,...,8 labels the eight eigenvalues given in (20,^) and those of the damping 
matrix T. Since we are interested in the correlation function of the operator A zz q, we have 
to calculate the element C\\ of the correlation matrix in the original basis, C = UCU^ . The 
transformation U being unitary, we have 

C n (z) =J2\Utj\ 2 C jj (z), (23) 

3 

where the vector of coefficients 

(\Uijl 2 ). = (0, 1/4, 1/4, 0, 1/8, 1/8, 1/8, 1/8) (24) 

satisfies the condition J2j \Uij\ 2 = 1- 

The damping rates Tjj are calculated to second-order in terms of the impurity-phonon cou- 
pling potential (||). Each rate is given as the convolution of an uncoupled correlation spectrum 
Cfi(uj) and the dissipation spectrum (|47|). The derivation is straightforward; we merely give 
the result 

r = dmg(ri,ri,r 2 ,r2,r3,r4,r3,r 4 ), (25) 

with 

r 1 = r (i + n(2A)) 

T 2 = T n(2A) 

r 3 = r (i + 4n(2A)) [Z0) 

r 4 = r (3 + 4n(2A)). 



Here we have used the Bose occupation factor 



™( 2A ) = ^n ( 2? ) 



and the coupled phonon Green function fl47|) evaluated at the frequency 2A/h, 

r = 2r»(2A/ R )=^/„||^. (28) 

The geometric factor f a arises from the different weight for longitudinal and transverse modes 
in the coupling matrix elements. 

Now we have an explicit expression for the frequency dependent function (|2^). After inverse 
Laplace transformation we obtain the time-dependent correlation function of the operator 

^>-zzCI> 

Cn(i) = - (mie- ri< + m 2 e~ r2< + m 3 cos(2At + 5 3 )e- r3t + m 4 cos(2At + <5 4 )e~ r4 *) , (29) 

which constitutes the basic theoretical result of this paper. 

As the most salient features we note the presence of two relaxation features with different 
relaxation rates Ti and r 2 . The remaining terms (the socalled resonant terms) show oscil- 
lations with twice the tunneling frequency A/h and different damping rates, with tan<5j = 
KTi/2A. The amplitudes m, vary significantly with temperature; in the high-temperature 
limit all of them tend towards unity, mi — ► 1, whereas for T — v one finds 7713 = 4 and 
m\ = m<i = ttt-4 = 0. Thus the relaxation contributions disappear at zero temperature. 

The spectra of the response function (|12j ) and the correlation function (|29| ) are related through 
the fluctuation-dissipation theorem, 

X "{uj) = (2/h) tanh(/37iw/2)C , ; , 1 (cu). (30) 

With the usual approximations for the hyperbolic tangent function we thus obtain 

X'V) = Phu [Li(0) + L 2 (0)] + tanh(/3A) [L 3 (2A) - L 3 (-2A) + L 4 (2A) - L 4 (-2A)] , (31) 

with the weighted Lorentzian 

^-TpcnlW (32) 

Clearly, the spectrum of the response function x"(w) contains the same physics as the corre- 
lation function Cn; hence the discussion below (|29|) applies equally well to (pl|). 

The tunneling impurities can be probed by an external ultrasonic wave with propagation 
vector k and polarisation vector e a . The attenuation of the latter is given by the internal 
friction 

7, 

whereas the variation of the sound velocity 



Qa^nh^^x"^), (33) 



5v a jl 
= nhk a - — 

v a 2gv_ 



— = nh ka ^ jX '(u;), (34) 



depends on the real part of the impurity dynamic susceptibility x' ■ Here n is the number 
density of the impurities; the geometric factor h^a, as defined in (|52|), accounts for the 
orientation of k and e a with respect to the crystal axes. 

Since the tunneling frequency A/% is of the order of 20 GHz, the second term of ( |3l|) is 
immaterial at acoustic frequencies; and the internal friction involves the relaxational part 
only, 



Qa L = nh ka >;: 



gvl k B T 



luTi UjTo 



mi „ 2 + m 2 2 2 



(35) 



Note that both amplitudes are exponentially small at low temperatures, m-i + m 2 = 2* 
sech(/3A/2) 2 , though m-i vanishes faster than m 2 . On the other hand, T± tends towards 
the constant Tq, whereas T 2 disappears. Thus in the intemediate range k-^T ~ A both 
contributions may be relevant. 

Because of the different temperature dependence of the two rates and the corresponding 
amplitudes, the variation of Q^ 1 with frequency changes very much with the ratio u/Tq. 
The origin of the terms in ( pop is quite obvious in view of the level scheme of Fig. 111. The 
contribution involving T 2 stems from relaxation between the states of the level Ti u ; the 
required intermediate transition to the top level A2 U accounts for the temperature factor of 
T 2 - On the other hand, the term with Ti results from the upper triple level T2 g ; relaxation 
occurs through the intermediate state Ai g ; phonon emission in this downscattering process 



gives rise to the temperature factor of T\. The rates of the resonant contributions in (31) 
can be discussed in a similar fashion. 

Finally we note that the relaxation rates T± and T2 satisfy the condition Ti = e@ 2 T 2 . They 
are related to the phase decoherence rates T 3 and T4 through Ti + T 2 = (l/4)(r 3 +T4); thus 
the transverse rates are larger than the longitudinal ones, contrary to the well-known ratio 
for the two-level system, \jT\ = 2/T 2 . 

3.2 Effects of asymmetry 

In order to obtain a realistic model for tunnel impurities, we have to account for the random 
asymmetry fields (p!0[). At sufficiently low concentrations, c < 100 ppm or n < 10 19 cm" 3 , 
typical values of Vi are significantly smaller than the tunnel energy A. Thus we may use a 
perturbation expansion in powers of the small parameter Vi/A. 



Each pole of the dynamic susceptibility (31) is characterized by an amplitude rrii, a damping 
rate Tj, and an oscillation frequency that takes the values ±2A or 0. It turns out that the 
lowest-order corrections to rrii and Tj are quadratic in Vi/A and thus are of little significance. 
Similarly, the finite poles at ±2 A are hardly affected by the small asymmtry energies Vi. 



A significant change occurs, however, in the zero-frequency poles. When evaluating (15) with 
both (||) and (|10| ) and calculating the lowest-order corrections in the diagonal representation 
(|22|) , we obtain the set of frequencies 



fl w diag(r], rj, -77, -rj, -2A, -2A, 2A, 2A) (36) 

with 



rt = 1/2 (vf-vi)/ A. (37) 



Accordingly, the relaxation contributions -^i(O) and £2(0) in (^Tj) have to be replaced by 
\(Li(rj) + Li(-rj)), with i = 1,2. 

Since experimentally frequencies are much smaller than the tunnel frequency A/h, the imag- 
inary part of the susceptibility is determined by the low-frequency poles ±77. When denoting 
the average with respect to the random fields V{ by a bar, we have 

PR = \f3u (T^jj + L^v) +W + Z2H7)) • (38) 

When inserting the inverse Fourier transform of the Lorentzians, the integral over V{ involves a 

Gaussian with complex width parameter (l/2<r 2 ) ±it/hA. Performing the Gaussian integrals 

we find 

wA f°° cos(wt) 



x "^ = f3 ^2j dt 7f^(, mie +m2e J' ( 39 ) 

with u) = ujA/cj 2 and T = TA/a 2 . We recall that we have used ksT, A 3> Twj, %T. 

Now we turn to the real part that is obtained from ( pOj) through the Kramers-Kronig re- 
lation. Proceeding as above we insert the inverse Fourier tranform; resorting to the usual 
approximations for the temperature-dependent factors and assuming f3rj/2 < 1, we obtain 
the susceptibility 

X'M = -^p- Z tanh (/3 A) + /> EE T ^^ + r ^ ^) 

Averaging over the distribution of asymmetries yields 

77 , 7713 + m4 1 , ,„.. v^ m; „ I ,, wA Z" 00 . sin(cDt)e _ri4 1 
x .(„) 1 «-L_l_ teDh()J A )+ g i ^|i__ r ^ *^)j^ (4!) 

3.3 Comparison with a two-level tunneling system 

When investigating the dynamic properties of substitutional impurities, various authors have 
used a two-state approximation. For this reason we briefly discuss the response function of 
a two-level system (TLS) and compare it with the present results for the eight-state system. 
With the notation introduced in Sect. 2, the Hamiltonian of a symmetric TLS reads as 

H = -Aa x + a z J2 7<*e§ + # P honon, ( 42 ) 

which is nothing but the one-dimensional version of @ and (||). When evaluating the dissi- 
pative two-state dynamics one finds that there is no relaxation contribution to the response 
function of a TLS; its susceptibility is of purely "resonant" character and similar to the second 
term of (|31~|), 

Thus the two-state approximation misses the observed relaxation peak in the dynamic re- 
sponse. In order to mend this failure, an asymmetry energy has been added to (|42|), yielding 



10 



H = 7}Aa x + ^ecr 2 , with a two-level splitting E = y/A 2 + e 2 . The spectrum of the dynamic 
susceptibility then reads as 

X 'L = * tanh( W 2) £ { ^^ %)2 + T2 + i-ch 2 (^/2)^ T M. (43 ) 

The amplitude of the relaxation feature is proportional to e 2 /E 2 and thus vanishes in the 
limit of zero asymmetry energy. This feature distinguishes the actual eight-state system from 
the two-state approximation. The former contains a strong relaxation feature, even at zero 
asymmetry energy. The two-state approximation, however, fails to account for relaxation in 
the symmetric case e = 0. 



4 Experimental detail and results 

The crystals measured here were seed pulled from the melt at the crystal growth facility of 
the Cornell Center for Materials Research. The starting powder of KC1 was Merck Industries 
"Superpure" grade with nominal impurity levels of less than 1 ppm. (A check of the OH 
level in our crystal of "pure" KC1, via UV absorption, confirmed a concentration of 0.5 ppm.) 
The cyanide added to the melt was taken from a seed-pulled single crystal of KCN for which 
the starting powder had been vacuum baked to remove H2O. 

The composition of the final crystal was analyzed via infra-red absorption (~ 2100 cm -1 ) to 
an estimated accuracy of 10% relative error, by comparing the area of the absorbance peak 
to a standard measured for us by Fritz Liity at the University of Utah where he measured 
the CN _ vibrational spectra at liquid nitrogen temperature and calculated the concentration 
based on the known oscillator strength of CN~ ]13| . 

The internal friction sample was cleaved directly from the IR absorption sample to ensure a 
known CN~ concentration. The internal friction was measured with a composite torsional 
oscillator as described in Ref. 14]. In this method, a 90 kHz quartz transducer (2.5 mm 



diameter) and the sample form a composite torsion bar. The quartz end is attached to a thin 



Be-Cu pedestal [14| by an approximately 25 mg drop of Stycast 2850FT epoxy, and since the 
KC1 crystals are quite fragile, a 0.25 mm indium foil was epoxied between the crystal and 
transducer to "cushion" the difference in thermal contraction rates of the two components 
upon being cooled. The sample length is tuned to be one half of a shear wavelength, so that 
the composite oscillator has a resonance frequency at room temperature within 1% of the bare 
quartz crystal resonance. This adjustment ensures that the epoxy and indium joint between 
the quartz and sample has almost zero strain, and therefore contributes minimally to the 
observed internal friction. (The epoxy and epoxy/indium junctions produce a background 
contribution to the internal friction of less than 10 -5 at low temperatures.) 

The oscillator is driven by a set of electrodes which form a quadrupole configuration around 
the transducer and which simultaneously drive and detect its motion. The internal friction 
of the sample Q~£ is determined from the quality factor Q comp of the composite oscillator 
resonance by 

'I sa + (l + a)I tr ' 
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Figure 2: The internal friction data of four KC1 samples with different concentrations of 
CN _ dopants are plotted versus the temperature. 



where It r and I sa are the moments-of-inertia of the transducer and sample respectively, and a 
is a correction for the attachment of the transducer to the thin Be-Cu pedestal, a ~ 0.06[14]. 
The relative change in speed of sound can be found from the change in resonance frequency 
of the composite oscillator 



Av 



V sa 



I sa + (1 + a)L, 



tr 



Isa 



f ' 



(45) 



where A/// is determined relative to an arbitrary reference frequency usually at the coldest 
temperature of the run. Details on evaluating internal friction from a torsinal oscillator can 
be found in Ref. fll5| ]. The measurements below 1.5 K were made in a dilution refrigerator, 
and those from 1.5 to 300 K in an insertable 4 He cryostat [f[q|. 

We have investigated a T2 9 -mode ( e = [110], k = [001] ) in three samples doped with 45, 370 
and 4700 ppm CN~, and a pure sample; the latter provides the background due to the host 
crystal and thus permits us to determine the impurity contribution to the elastic response 
function. 

Figure |2| shows the internal friction as a function of temperature, cyanide concentration 
ranging from to 4700 ppm. The data consist of several peaks. The dominant feature 
occurs at about 450 mK in the 45ppm-data and moves to about 350 mK in the 370ppm- 
data; in the strongly doped sample (4700ppm) it results in a broad shoulder beyond the 
lowest temperatures investigated. At higher T there are additional peaks with much weaker 
intensity, which may well result from strongly coupled defect pairs |18fl . The internal friction 
of the undoped sample is significantly smaller, although still remarkably high. (Even after a 
check for purity, there is still some unexpected damping in the "pure" KC1. The expected 
background internal friction value for pure crystals using this technique is less than 3 x 10~ 6 
below 1 K, as measured in quartz |15|.) 
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Figure 3: Relative change of sound velocity Sv/v for four KC1 samples with different con- 
centrations of CN~ dopants are plotted versus the temperature. For convenience, the data 
are multiplied by the factors of 10. Note the vertical shift by steps of 0.5 with respect to the 
somewhat arbitrary origin at zero temperature. 

The relative change of the sound velocity is plotted in figure ||. The temperature dependence 
of the data may be decomposed in two features. First, there is a contribution proportional 
to tanh(E/2kT) with an energy E of a few K, that is characteristic for quantum defects 
and shows the well-known 1/T-dependence at higher temperature. Second, the minimum at 
T = 0.7 K indicates a relaxation process. 

Though the minimum broadens with increasing doping, the three samples show quite a similar 
behavior, and the change in sound velocity is roughly linear in the impurity concentration. 



5 Comparison of theory and experiment 

We focus on the data of KCL with 45 ppm CN~ since our theory is correct for low doping 
only, i.e. for impurities without interaction. With increasing doping, however, their dipolar 
interactions are no longer small. The collective dynamics results in the broad internal friction 
spectrum shown in figure || at higher concentrations; the sharp features indicate strongly 
coupled impurity pairs. 

The theoretical expressions for symmetric tunnel impurities involve the tunnel energy A and 
the damping rate Tq (p8[); taking also finite asymmetry energies as in Sect. 2.3 into account, 
the Gaussian width a provides one more parameter. Both the rate Tq and the prefactors 
of the internal friction and the change of sound velocity (|3"3| |34|) depend on the deformation 



potentials j a , the sound velocities v a , and the mass density of the host crystal g. The 
latter quantities are well known; since the longitudinal sound velocity is significantly larger 
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Figure 4: The [lll]-defect theory is fitted at the data of the internal friction of the KC1: 
45ppm CN~ sample. The full line is a fit including asymmetry effects and the short-dashed 
line without these. The long-dashed line shows a plot of the internal friction of a two-level- 
system where we set (e/E) = 1 whereas all other parameters was choosen the same as in the 
former fit. 



than the transverse one, we have v^~ 3> vT and retain transverse sound waves only, with 
v t = 1.7 km/s and g = 1.989 g/cm 3 . (We neglect the weak dependence of the sound velocity 
on the propagation direction.) 

The tunnel energy A/ks of CN in KC1 has been measured by various techniques, with differ- 
ent experimental results; paraelectric resonance: 1.87 K |2^| ; excitation of optical vibrations: 
1.73 K |p3[ |; specific heat: 1.6 K f24[| . The relatively large discrepancy between these values 
may be due to experimental uncertainties and different parameters of the samples, in partic- 
ular impurity concentration. From our fits to the 45 ppm data we obtained a value of 1.55 K, 
which is in reasonable agreement with the above results. Regarding the asymmetry energy, 
we use a/ks = 0.0145 K for the width of the Gaussian distribution. Finally, the transverse 
deformation potential jt constitutes the most important parameter, appearing both in the 
prefactors and relaxation rates. From our fits we obtain the value jt = 0.192 eV which is of 
the same order of magnitude as the values obtained for OH impurities in KC1 (1 eV) |^5| , in 
NaCl (1.0 eV) @ and OD in NaCl (0.34 eV) ||. 

First we discuss the internal friction as shown in figure [|. The triangles are the background 
corrected data of the 45 ppm sample, i.e., we have subtracted from the original data the 
values for the undoped sample. The full line is a fit including asymmetry effects ([^) where 
we used the parameters given in Table 1. 

The temperature where the maximum occurs depends strongly on the tunneling amplitude 
A and the fraction uj/Fq; the maximum value of Q _1 varies with the prefactor and the 
asymmetry width as far as the mean splitting is bigger than the frequency f\ = a 2 /A > uj. 
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Figure 5: The [lll]-defect theory is fitted at the data of the relative change in speed of sound 
of the KC1: 45ppm CN~ sample. The full line is a fit including asymmetry effects and the 
dashed line without these. 



The relaxation peak arises where the external frequency is comparable to the slow rate T2, 



and T max is determined by the relation T2(T„ 



10. Thus the exponential decrease of T2 



is essential for the existence of the relaxation maximum. 

Note that this mechanism does not exist for two-level tunneling systems, even with large 
asymmetries, since their relaxation rate tends towards a constant Tq. As a consequence, for 
U) <C To the relaxation maximum is suppressed by a factor uo /Tq, and the internal friction 
of a two-level system is by two orders of magnitude smaller than that obtained for the 
[111] -impurity model. Note that the parameters gathered in Table 1 result in a rate constant 
To = 1.35 x 10 9 sec -1 which is indeed much larger than the frequency omega = 5.3 x 10 5 sec -1 . 
For systems with 77 > u> the peak is suppressed by oj/tj (see eq,(|38|)) which means that the 
peak contribution stems from systems with small asymmetries. Therefore the peak in the 
internal friction data determines only a combination of the prefactor and the asymmetry 
width. 

Our simple model accounts well for position and height of the peak and for its shape close 
to the maximum. Yet it strongly underestimates the wings at temperatures well above and 
below T max . As to the excess spectral weight observed at low temperatures, a more realistic 
distribution of asymmetry energies would probably give a better agreement than the Gaussian 
used in (|39l). The additional peak at about 2 K may well be due to strongly coupled pairs. 

Now we turn to the relative change of sound velocity as shown in figure [|. The triangles are 
data for the 45 ppm sample. The dashed line is a fit without asymmetry, given by (p4|), and 
the full line with asymmetry (|4l|). The asymmetry distribution with width a merely smears 
out the temperature dependence. Thus the fit of Fig. [5] mainly depends on the tunneling 
amplitude and the prefactor but it hardly varies with the relaxation rate and the asymmetry 



15 



uj/2tt (sec l ) 


n (cm 3 ) 


q (gem 3 ) 


^i (km/s) 


A/fc B (K) 


It (eV) 


a/fe B (K) 


84352 


1.7 x 10 17 


1.989 


1.7 


1.55 


0.192 


0.0145 



Table 1: Parameters used for the fits of Figs. || and[B| 

width. 

Regarding both Figs. || and [B|, our model fits remarkably well the absolute values of Sv/v 
and Q _1 , the hump in the sound velocity, and the position of the relaxation peak in the 
internal friction. The transverse deformation potential is the most relevant parameter that 
determines both the absolute values of Sv/v and Q _1 , and the relaxation maximum of the 
latter; we find the value 7j = 0.192 eV. (Note the factor 2 in our definition of the phonon 
coupling (||).) 

The width of the relaxation peak of Q~ l and the low-temperature wing indicate the relevance 
of the asymmetry energies (IToT) . The poor quality of our fit at low T may well be due to the 
choice of a Gaussian distribution for the asymmetries. Nonetheless, our results confirm that 
for the 45 ppm sample, the mean asymmetry is by two orders of magnitude smaller than the 
tunnel energy, which is consistent with the model assumptions. 

The present theory accounts fairly well for available data on the relaxation behavior of [111]- 
impurities in alkali halides investigated with T2 9 -modes, in particular with respect to the 
temperature dependence. Yet the multiexpenontial decay with the two rates T\ and T2 
gives rise to quite an intricate relaxation spectrum; thus an experimental investigation of 
the frequency dependence would seem most promising. A simliar behaviour was also found 
for an E g -mode || contradicting the simplest model of an [111] -defect. Even asymmetries 
can not explain this feature in a simple way. Another open question concerns the thermal 
conductivity Pj. The relaxation feature obtained in the present work should significantly 
contribute to thermal resistivity at low T. Still, the spectrum of the defects obtained from 
the data |2| suggests a more complicated elastic response spectrum and may well require large 
asymmetry energies. 

Another interesting point would be to extend the present work to the case of strong doping; 
this is certainly not an easy taks, given the serious difficulties encountered already in the 
treatment of interacting two-state systems (see, e.g., Ref. [pTJ] .) Note, however, that our sound 
velocity data in Fig. || depend, roughly speaking, linearly on the impurity concentration. 
Thus it would seem that they do not fullfil the strong-coupling criterion of ||, i.e., the 
dipolar interactions have not yet destroyed the coherent tunnel motion. 

6 Summary 

We have investigated the relaxation of an impurity ion in alkali halides arising from the 
coupling to elastic waves. We briefly summarize the main results. 

(i) The various elastic and inelastic phonon-mediated transition between the eight quan- 
tum states give rise to an intricate temperature and frequency dependence of the relaxation 
contributions to the internal friction and the sound velocity (see Fig. ffl.) Unlike two-state 
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tunneling systems, [lll]-impurities show two relaxation rates. At low T, the smaller rate 
decreases exponentially, T2 = e~ 2A ' kT ro, whereas the larger one tends towards a constant, 
Ti = r . This corresponds to a multiexponential decay of the time-dependent response and 
correlation functions (O) and (p9|). 

(ii) A most particular relaxation behavior arises for external frequencies u> that are smaller 
than the constant rate, uj < Tq. Then the relaxation maximum occurs at a temperature 
where the smaller rate 1?2 is equal to the external frequency u. This maximum is the more 
relevant as at low T the spectral weight of the slow contribution exceeds by far that of the 
faster rate Ti. 

(iii) Comparison with recent data on KC1:CN proves the relevance of this relaxation mech- 
anism. The exponentially decreasing rate T2 explains the large amplitude of the relaxation 
peak shown by the external friction Q~ l as a function of temperature at v = 84 kHz. By the 
same token, our model provides a good description for the hump in the sound velocity. 

(iv) The prefactors of 5v/v and Q^ 1 and the relaxation rates are related by the values for 
the deformation potential 7, sound velocity v, and mass density p. By taking 7 as a free 
parameter, we obtain for reasonable values of 7 (compared with simliar materials) satisfying 
fits for 5v/v and Q _1 , involving both absolute values and the temperature dependencies. 

(v) The phenomena mentioned in (iii) cannot be explained in terms of relaxation of corre- 
sponding two-level systems. For the latter, the low-temperature rate tends towards a constant 
and thus may exceed uj at any T, whereas the small rate T2 of a [lll]-impurity inevitably 
meets uj at some T and gives rise to a relaxation peak. Though certain aspects of the ther- 
mal and dielectric properties of such doped crystals are described by a ensemble of two-state 
systems, such a model fails in view of the acoustic properties, due to the multiexponential 
decay of the elastic response function ( |l~2|) . 
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A Geometric factor of the ESS 

Time evolution of the phonon heat bath is given by the lattice response function (t > 0) 

r(t) = E if <[<§(*),*§]>, ( 4 6) 

where i 7^ j and the bar indicates the average over crystal axes. The entries of the damping 
matrix ( |l8| ) are determined by the coupled phonon spectrum; in the limit of long wavelengths 
one has 

T» = ^ E ft^a-^—iS^ - ^a) ~ 5(U, + U, ka )], (47) 
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where 

f^ = (efi^k + e^k)k) 2 (48) 

accounts for the orientations of wave and polarizations vectors with respect to the crystal 
axes. For an isotropic phonon density of states, the sum over phonon modes becomes 

£/te (•••) — ^ 



£ £/£?(-) — i£H d ^ 2 (-) ( 49 ) 



where we have for each polarization defined the average value 

In the Debye model with Wk a = v a k, the damping spectrum finally reads as 

Since an elastic wave couples via each of the quadrupole operators to the defect according 
to (||) , the geometric factor hk a of the internal friction and the relative change of the sound 
velocity of an elastic wave with propagation vector k and polarisation vector e Q gets 

fc» = /£° + /£? + /£> (52) 

The factors f a are easily evaluated after expressing the unit vectors k = k/fc and e in polar 
coordinates 9 and (j). Putting k = (sin#sin</>, sin#cos</>, cos#), one finds for the longitudinal 
case e 1 = k the factor fi = A/15. 

Regarding the transverse modes, we obtain f t2 = 8/45 for e 2 = (cos <j>, — sin (j), 0) and ft 3 = 
2/9 for the remaining polarization e 3 . Though these values depend on the choice of e 2 and 
e 3 , their average ft = (1/2) (/j 2 + ft 3 ) = 1/5 is independent of the basis. Defining moreover 
the mean value of the three polarisation directions, / = (1/3) (fi + 2ft), we finally have 

4 12 

fi = —, ft = ~, / = -• (53) 

n lg , jt 5 , j g v ) 

Thus the quantity ( f48|) for longitudinal modes is by a factor fi/ft = 4/3 larger than for the 
transverse ones. In the average, propagation and polarisation directions are uncorrelated, 
resulting in 

/=((4 2 )(S 2 ) + (^ 2 )(e?))=2(l/3) 2 (54) 
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